Anti-correlation between energy-gap and phonon energy for cuprate Bi2212 

superconductor 



o 
o 

(N 

as: 

< 

o 

m 



c 

o 
o 

I 

5-1 

Oh! 



c 

o 
o 



> 
m 

O 

o 
o 



X 



Wei Fan 

A"e?/ Laboratory of Materials Physics, Institute of Solid State Physics, 
Chinese Academy of Sciences, 230031-Hefei, People's Republic of China 
(Dated: February 1, 2008) 

Using electron-phonon mechanism, we explains very well the spatial anti-correlation between the 
energy-gap and the energy of phonon mode for cuprate superconductor found in d 2 I/dV 2 spectrum 
by STM measurements [Jinho Lee, et.al Nature 442,546 (2006)], which is is the direct effect of a 
relationship M{uj 2 )\ = const. We calculate Tc maps on A — Qp plane which helps us understanding 
the relation M{ui 2 )\ = const and the superconductivity of superconductors. 

PACS numbers: 74.20.Fg, 74.25. Kc, 74.62.-c, 74.72.-h, 74.81.-g 



Scanning Tunneling Microscopy/Spectroscopy 
(STM/S) has been used to explore surface struc- 
tures and electronic structures of superconductors 
The electronic inhomogeneity such as the inhomoge- 
neously distributed energy-gaps and local densities of 
states has been found in High Temperature Supercon- 
ductor (HTSC) @,[3. The intrinsic inhomogeneity even 
exists when translational invariance of lattice structure 
preserves. If total momentum of cooper pair isn't equal 
zero, in real-coordinate space the distribution of cooper 
pairs is inhomogeneous. The Fulde-Ferrell-Larkin- 
Ovchinnikov (FFLO) phase can be realized by applying 
magnetic field or stress field to superconductors. 

An interesting property of electronic inhomogeneity of 
superconductor by the observations of STM/S for HTSC 
Bi2212 is the spatial anti-correlations of energy-gap and 
energy of boson mode [3j. The understanding of origin 
of the anti-correlation is very useful to understand the 
microscopic mechanism of superconductivity for HTSC 
materials. The boson mode for Bi2212 has been identified 
as phonon mode due to the isotope effects. 

The out-of-plane buckling modes, in-plane breath- 
ing modes and the phonon modes related apical oxy- 
gen atoms are expected to contribute to superconduc- 
tivity of HTSC materials. The out-of-plane buckling 
modes coupling with the electrons near anti-node region 
have small momentum transfer and significant band re- 
normalization effects near transition temperatures which 
had been found in ARPES spectrum and theoretical cal- 
culation 

In this paper, we preform the standard calculations 
of Nambu-Eliashberg strong coupling theory and focus 
on an important relation M(u> 2 )\ = r\ = const. The 
superconductor parameter r\ = p(0)(J 2 ), the product 
of the electronic properties p(0) the density of state at 
Fermi energy and the electron-phonon matrix elements 
J, characterizes the chemical environments of atoms and 
almost keeps as a constant against the simple structural 
changes such as the isotope substitution. We assume 
the electron-phonon mechanism is pairing mechanism for 
superconductor Bi2212 in this paper because of the iso- 
tope effects and identification of phonon mode in ARPES 
spectrum [f|, STM/S experiments Q. In terms of the 



relation M(uj 2 )X — n, we explain the anti-correlation be- 
tween energy-gap A and phonon energy £1 p . We also plot 
two Tc maps of transition temperature on A — ilp plane, 
which are very useful to understand superconductivity of 
conventional and unconventional superconductors. 

The theoretical methods used in this paper are (1) 
the real-energy method by having analytically contin- 
ued the imaginary Eliashberg Equation to real axis 0, Q 
and (2) the imaginary-energy Matsubara method [9( of 
strong coupling superconducting theory under isotropic 
approximation. After having resolved isotropic energy- 
gap equation, we can obtained the anisotropic energy- 
gap by directly multiplying isotropic energy-gap A(u) 
with the anisotropic function such as for the d-wave the 
anisotropic energy-gap A(w)cos(2#) and then averaging 
it over 9 angle. The anisotropy of energy-gap comes from 
the anisotropic electron-phonon interaction [3, Q . The 
real-energy method gives the results of densities of states, 
energy gaps and Tc- The advantage of imaginary meth- 
ods is that we can quickly obtain the T^ compared with 
the real-energy methods. 

The densities of state of phonon is calculated by fol- 
lowing formula 0] 



F(u) 







others, 



(1) 



where f2 p is the energy of phonon mode, u; 2 the half-width 
of peak of phonon mode and ui^ = 2ui 2 . The effective 
electron-phonon interaction a"^{v) is defined by 

a 2 y (u,)F y (u,) = - / V - \ ^B y (p-p',w)\Dl_ p , >y \\{2 

W J s V F Js> V F> 

where W = l/(27r) d_1 J s d d ~ 1 p/vp and v F Fermi ve- 
locity. In 2-dimensional limit the coupling between dif- 
ferent layers being ignored, we can solve the Eliashberg 
equation for 2-dimensional systems. Compared with 3- 
dimensional system, the energy-gap equation keeps un- 
changed and all changes arc included in the Eliashberg 
function a 2 F(u>). Complete treatments of anisotropic ef- 
fects need the k-dependent energy-gap function A(fc, u)) 
and re-normalized factors Z(k,ui). 
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FIG. 1: (a). The anti-correlation between energy-gap A and 
the phonon energy tip and corresponding isotope effects for 
16 substituted by 18 0. (b,c). The s-wave and d-wave densi- 
ties of states at 5 K and 80 K. The definition of energy gap 
A by densities of states (DOS) and energy gap with d-wave 
symmetry in (b). 



We can write the parameter of electron-phonon 
interaction A = ^2 y \, where parameter X y 
1 duJctyFy(uj)/uj is the contribution of y mode to 
total A. The moments (uj n ) of the distribution 
function (2/ X)^2 y 2a 2 F y (uj)/cj are defined as (lu 71 ) = 

2/A J2 y Jo°° dujayFy^uj™^ 1 . In this work we assume 
that a y is constant around each peak in phonon 
spectrum. Under this approximation, we only con- 
sider a 2 F(oj) as a whole, the so-called Eliashberg 
function with above Lorentz shape. In practice, by 
modifying the height of peak of Eliashberg function 
a 2 F(uj), we can control the change of the param- 
eter A of electron-phonon interaction. The densi- 
ties of states of electron are obtained from the for- 
mula p(uj) = pN{(jj)Re(uj I \Joj 2 — A(w) 2 ). If the 
anisotropy of energy gap is considered such as with 
the d-wave symmetry, the density of state p(uj) = 
PN^iReiuj/^/uj 2 - (A(w)cos(0)) 2 )) e , with the (...}„ 
describing the average over 9. 

We choose the phonon spectrum of Eq. {1} with 
phonon-energy about Op=52 meV, which is the mean 
phonon-energy observed in the tunneling spectrum of 
Bi2212 0. The P distributes from 30 meV to 70 meV 
and within the range from about 36 meV the Bi ff phonon 
coupling to electrons at antinode region and to 70 meV 
the in-plane breathing mode coupling with electrons at 
node region based on the observation of ARPES spec- 
trum [6j. 

If the parameter A of electron-phonon interaction is 2.0 
and the Coulomb pseudo-potential p*—0.24, the transi- 
tion temperature Tc is about 90K close the experimen- 
tal values. From STM experiments of Bi2212, the en- 
ergies of phonon have a distribution with width about 
2^2=16 meV. The energy gap A and the A the parameter 
of electron-phonon interaction should have distributed 
with their widthes respectively. If we consider a material 



with multi-species atoms, total electron-phonon interac- 
tion parameter should be expressed as A = ^ TO A m with 
A m = Vm/Mm(^rn) the electron-phonon interaction pa- 
rameters for every species atoms. If the phonon modes 
correlated with the vibrations of oxygen atoms have most 
significant contributions to superconductivity, we can get 
simple relation M(co 2 )X = Tj = const. Additionally, we 
assume it is correct for isotope substitution, when mass of 
oxygen atom increases from 16 a.u to 18 a.u. Both A and 
(ijj 2 ) should be changed to keep M(ui 2 )X = r\ = const. 

The Fig. fl} plots the A — O relation for the sam- 
ples with 16 O and 18 O. Certainly, the products M{uj 2 )\ 
are all the same for all the points in the Fig. (JTJ (a) . 
The curve for 18 are beneath the curve for 16 0. The 
very similar result is found in Fig. 5(a) in Ref. 
The energy-gaps decrease with increasing phonon en- 
ergy Op because the parameters of electron-phonon in- 
teraction increase with decreasing phonon energy because 
M{ui 2 )X = const. In a sample of superconductor, the in- 
homogeneously distributed phonon-energies Op(r) lead 
to the in-homogeneously distributed energy-gap A(u>, r), 
where Op is large, A(w,r) small. The vertical shift of 
A — Op curve downward to smaller value for larger mass 
is due to the increasing parameter of electron-phonon in- 
teraction with decreasing Mass because of the relation 
M(uj 2 )X = const. However, the energy gaps A in this 
work are only about half of values in Ref. Q and the verti- 
cal shift of 4.0 meV in A smaller than 5.6 meV in Ref. [![. 
The reason is that the energy gaps in this paper are 
directly correlated to the formation of coherent cooper 
pairs and the energy gaps in Ref. [|[ are the mixing of 
superconducting gap and pseudo-gap of spin fluctuations. 
The values of superconducting energy gap is about from 
23 meV larger than the values 19-20 meV in this work. 
The effects of the relationship M(uj 2 )X = const had al- 
ready been discussed in Ref. [{J, [l(| ■ 

A Tc map on the parameters space A — O — u* is 
very helpful to understand the superconductivity with 
electron-phonon mechanism. The important role of the 
relationship M(lo 2 )X = const can be shown clearly on 
the Tc map. In this paper we provide Tc maps on 2- 
dimensional A — O plane with 50 x 50 mesh-grid with 
parameters A from 0.5 to 4.0 and Op from meV to 80 
meV. We plot two different Tc maps with the Coulomb 
pseudo-potential /i*=Q.l Fig. © and 0.24 Fig. ©. For 
every mesh points we can calculate a Tc using imaginary- 
energy Matsubara methods. An important point is that, 
for high energy phonon mode which is the vibrations of 
relatively light oxygen atoms, only intermediate A(< 2.0) 
or intermediate electron-phonon coupling interaction can 
guarantee the values of Tc within the range of HTSC 
materials. As we expect that for large Coulomb pseudo- 
potential p* =0.24 the larger parameters A of electron- 
phonon interaction are needed to obtained the same Tc 
for the general values p*=0.1. 

In two figures we also plot two curves of OpA=320 
and 3600 (meV) 2 . We have been ignored the changes 
of mass of atoms in the relation M(ui 2 )X = const. We 
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FIG. 2: The Tc map on A - li p plane with fx* =0.10 ob- 
tained from imaginary-energy Matsubara method. The bold 
solid line and dash-dot line are the curves Af2p=320 and 3600 
(meV) 2 . The region A includes the contour lines from K to 
30K, B from 30 K to 50 K, C from 50 K to 100 K and D from 
100 K to 150 K. 
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FIG. 3: The Tc map on A - O p plane with //=0.24 ob- 
tained from imaginary-energy Matsubara method. The bold 
solid line and dash-dot line are the curves Af2p=320 and 3600 
(meV) 2 

can see that the two curves cross smaller number con- 
tours lines when A > 2.0 with small phonon energies. 
This fact indicates that product {w 2 )X is an important 
parameters characterized Tc. Additionally, the curve for 
f2pA= 3600 (meV) 2 enters the region with higher tran- 
sition temperature compared with the curve for flpX— 
320 (meV) 2 . Although the range of A > 2.0 is unfavor- 
able superconductivity due to the lattice instability, the 
product is still good parameters characterized Tc . When 
superconductor under small structural modifications or 
isotope substitution, the corresponding parameters A and 
Hp move alone these curves because the chemical envi- 
ronments keep almost unchanged. 

For high-energy phonon, the smaller change of A can 
induce large change of Tc because the curve fJpA=3600 



cross more contour-lines than for low-energy phonon. 

TABLE I: Transition temperatures Tc calculated from the 
real-energy method T^ e and imaginary-energy Matsubara 
method Tg™. The relative errors e = {T§ e -T c m ) /Tg° , where 
Tc" = (Tg e + T c m )/2 





Q p (meV) 




A 




T c m (K) 


e 


(1) 


25 


0.10 


0.8 


14 


15 


0.069 


(2) 


25 


0.10 


1.5 


33 


37 


0.114 


(3) 


25 


0.24 


1.8 


34 


39 


0.137 


(4) 


25 


0.24 


2.0 


38 


43 


0.123 


(5) 


52 


0.10 


0.8 


33 


35 


0.059 


(6) 


52 


0.10 


1.5 


79 


82 


0.037 


(7) 


52 


0.24 


1.8 


79 


86 


0.085 


(8) 


52 


0.24 


2.0 


90 


95 


0.054 



So the superconductivity correlated with high-energy 
phonon is structural sensitivity. The energy of phonon 
plays more important role for high-temperature super- 
conductor. Because the atoms with smaller mass have 
larger frequency of vibrations, almost all superconduct- 
ing materials with high transition temperature contain 
light atoms such as Carbon atoms in fullerides and oxy- 
gen atom in HTSC materials. However these materials 
face the problem of instability of superconductivity in- 
duced by the structural instability. 

The Tc values obtained using two different methods 
(real-energy and imaginary methods) only have small dif- 
ferences with relative errors e < 0.14 if the high-energy 
cutoff is larger than 1000 meV and the low-energy cut- 
off is smaller than 1.20 meV. In the real-energy method, 
the energy cutoff is about 1000 meV at high energy and 
1.20 meV at low energy. In imaginary method, we use 
about 200 Matsubara energies to resolve the Eliashbeg 
equation. From the table (|TJ the larger relative errors 
happen at strong coupling regime with larger Coulomb 
pseudo-potentials and however with smaller phonon ener- 
gies such as the case(2,3,4). With increasing the energies 
of phonon from 25 meV to 52 meV, the relative errors 
will smaller than 0.1. The absolute values of errors from 
5K to 7K for the case (6,7,8) are small for high tran- 
sition temperature Tc however large for the case (3,4) 
with relative low transition temperatures. From above 
comparisons, we can say that the Tc map in Fig. @ 
is helpfully although there are some uncertain from the 
numerical calculations. 

In summary, we find the anti-correlation between 
energy-gap A and phonon energy Qp is deduced from 
the simple relation M{w 2 )\ — const in strong coupling 
theory. We also plot two Tc maps in A — Qp to under- 
stand the above relation and the superconductivity. This 
work is supported by Knowledge Innovation Program of 
Hefei Institutes of Physical Sciences, Chinese Academy 
of Sciences, and in part by National Science Foundation 
of China for CAS projects. 
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